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We reveal a close relationship between quantum metrology and the Deutsch-Jozsa algorithm on continuous vari- 
able quantum systems. We develop a general procedure, characterized by two parameters, that unifies parameter 
estimation and the Deutsch-Jozsa algorithm. Depending on which parameter we keep constant, the procedure 
implements either the parameter estimation protocol or the Deutsch-Jozsa algorithm. The parameter estimation 
part of the procedure attains the Heisenberg limit and is therefore optimal. Due to the use of approximate nor- 
malizable continuous variable eigenstates the Deutsch-Jozsa algorithm is probabilistic. The procedure estimates 
a value of an unknown parameter and solves the Deutsch-Jozsa problem without the use of any entanglement. 

PACS numbers: 03.67.Ac, 42.50.Ex, 06.20.-f 



I. INTRODUCTION 

Quantum metrology promises many advances in science and 
technology. Continuous variables (CV), i.e. eigenstates of an 
operator x with continuous spectrum, are natural candidates 
for optical implementations of quantum metrology protocols 
[1-3]. The importance of continuous variables for quantum 
metrology stems from the unconditional and efficient charac- 
ter of CV preparation, manipulation and detection techniques 
[4, 5]. In this paper, we devise an optimal parameter estima- 
tion procedure for continuous variables. Our procedure em- 
ploys a single continuous variable and estimates a value of an 
unknown parameter with Heisenberg-limited precision. Fur- 
thermore, for a particular, fixed value of the parameter in ques- 
tion the procedure behaves as the Deutsch-Jozsa algorithm for 
CVs. In fact our protocol extends the Deutsch-Jozsa algo- 
rithm over continuous variables presented by Pati and Braun- 
stein [6] . Instead of idealized, nonnormalizable and unphy si- 
cal states, we employ Gaussian states to represent continuous 
variables. Moreover, we define Gaussian states on a finite do- 
main, thus removing an unphysical, infinite speed-up over any 
classical procedure offered by the idealized states. An exten- 
sive analysis of the Deutsch-Jozsa algorithm over continuous 
variables was given by Adcock, H0yer, and Sanders [7]. 

The Deutsch-Jozsa algorithm is one of the first quantum 
algorithms, preceded only by the original Deutsch algorithm 
[8]. Even though the Deutsch-Jozsa problem is rather artifi- 
cial, the algorithm drew enormous attention due to the compu- 
tational speed-up over any classical procedure. The structure 
of the algorithm is simple enough to determine the source of 
this speed-up. The quantum superposition principle and con- 
sequent quantum parallelism that lie at the heart of quantum 
mechanics allows for the interference of many distinct com- 
putational paths, and allows the correct answer to the problem 
emerge in a single query. In other words, the Deutsch-Jozsa 
algorithm probes a global property of an unknown function 
f(x) and returns the result in a single run. 
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The paper is organized as follows. In Sec. II, we review 
the Deutsch-Jozsa algorithm over continuous variables and 
present its simplified version. In Sec. Ill, we review the ba- 
sic concepts in quantum metrology. In Sec. IV, we introduce 
a general procedure that unifies parameter estimation with the 
Deutsch-Jozsa algorithm, and analyze it in detail. Finally, we 
make some concluding remarks in Sec. V. 



II. DEUTSCH-JOZSA ALGORITHM OVER CONTINUOUS 
VARIABLES 

The generalization of the Deutsch-Jozsa algorithm to continu- 
ous variables was devised by Pati and Braunstein [6]. This 
generalization was implemented with idealized continuous 
variables defined on an infinite domain. However, we need 
to stress that any practical continuous -variable implementa- 
tion of the Deutsch-Jozsa problem can be realized only on a 
finite domain. Nevertheless, for simplicity, we first recall the 
Deutsch-Jozsa algorithm over continuous variables as origi- 
nally stated in Ref . [6] . 

The object of the Deutsch-Jozsa problem is to determine 
whether some function f(x) is constant or balanced. This is 
achieved by Alice and Bob playing the following game. Al- 
ice submits a value of x from — oo to +°o to Bob. Then Bob 
evaluates f(x), which can take only two values: or 1. Bob 
also promises Alice to use either balanced or constant func- 
tions. A constant function is either or 1 for all values of 
x G (— °°, +°°). A balanced function is for half of the values 
of x, and 1 for the remaining values of x. This is defined in 
terms of the Lebesgue measure jl on R: ji(x G M|/(jc) = 0) = 
fi(x G M|/(jc) = 1) [6]. The goal of this game is the same as 
the objective of the Deutsch-Jozsa problem, i.e. to establish 
if the function used by Bob was constant or balanced. Clas- 
sically, Alice would have to submit infinitely many values of 
x to learn the global property of f(x) with certainty. How- 
ever, if Bob can use a unitary black-box operation to calculate 
function f(x), then only a single function evaluation reveals 
the global property of f(x). In the setting of idealized contin- 
uous variables this would imply an infinite speed-up over any 
classical procedure. 

Let us now introduce in some detail the ideal Deutsch-Jozsa 
algorithm over continuous variables shown in Fig. 1. This 
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FIG. 1 : A quantum circuit representing the Deutsch-Jozsa algorithm 
over continuous variables. The quantum network Ndj consists of the 
Fourier transforms F and controlled black-box gate Uf applied to 
the register and target CVs prepared in the idealized position eigen- 
states |jto) and |tt/2), respectively. The last operation is an inverse 
Fourier transform that enables the interference of different computa- 
tional paths. 

implementation of the Deutsch-Jozsa algorithm employs two 
continuous variables, the so-called register and target CVs. 
Alice stores her query in the register CV and the target CV is 
used by Bob during function evaluation. The register CV is 
prepared in the position eigenstate |xo) and the target in the 
eigenstate |tt/2). The quantum network Ndj implementing 
the Deutsch-Jozsa algorithm is given by the following unitary 
transformation 

N D j = F~ l UfF r F t , (1) 

where F is the Fourier transform and r, t indicates the register 
and target CV, respectively. The Fourier transform applied to 
a CV in some position eigenstate \x) creates a superposition 
of all position eigenstates according to 

1 f°° 

F\x) = -y= J_Jye 2ixy \y), (2) 

where we used photon number units in which ft = |. The 
unitary black-box operator Uf evaluates a value of function 
f(x) and stores it in the state of the target CV: \x) \y) — » |x) \y + 
f(x)). Let us analyze the Deutsch-Jozsa algorithm step by 
step: (i) prepare the register and target CVs in an ideal position 
eigenstate |jco) and |tt/2), respectively; (ii) apply the Fourier 
transform F to the register and target CVs 

\s) =F r F t \x )\n/2) = - f°° dxdye 2ixx ^ iK y\x)\y)- 

71 J-oo 

(iii) following the action of a unitary black-box operator Uf 
the state of the CVs is given by 

1 f'°° 

u f\ s ) = ^J dxe 2ixx »e- iKf ^\x)F t \7l/2)- 

(iv) the quantum network Ndj is finalized with an inverse 
Fourier transform F~ l applied to the register CV. Therefore, 
the state of the CVs can be written as 

F~ l Uf\s) = - T dxdx' e 2ix ^-^e- in ^ x \x')F t \n/2)\ 

71 J — oo 

(v) following the quantum network Ndj, the property of func- 
tion f(x) is determined by projecting the state of the register 
CV onto the original position eigenstate |xo). The continuous- 
variable projection operator for idealized states can be written 

as 

rx +£ 

Px = dy\y)(y\, (3) 

Jxq — £ 



where £ is the spread around xo value, i.e. the continuous- 
variable measurement cannot be performed with infinite pre- 
cision. The orthogonal complement of P XQ is given by 

rx +£ 

P* =I-P X0 =I- dy\y)(y\. (4) 

Jxq — £ 

By construction, a complete set of orthogonal projectors P m 
satisfy the completeness relation Y^m^m = I and P m P m i = 
S mm fP m . If f(x) is constant then the measurement statistics 
based on the above set of orthogonal projection operators and 
assuming £ — >> is given by 

p(xo) = Tr[P XoP Dj] = 1, (5) 
P (xq) = Tv[P XoPd j} = 0, (6) 

where p(xo) is the probability of measurement outcome to be 
xq, p(xq) is the probability of a measurement outcome differ- 
ent than xo and pDj = Ndj \r) \t) (t \ (r\N^j . Conversely, if f(x) 
is balanced then the measurement statistics assuming £ — » is 
given by 

p(xo) = Tr[P Xo pD/]=0, (7) 
p(x ) = Tt[P XqP dj] = 1. (8) 

Therefore, if the state of the register CV remains unchanged 
then the function f(x) is definitely constant, and if the state 
of the register CV is not |jco) then the function f(x) is bal- 
anced. A single function evaluation solves the Deutsch-Jozsa 
problem. 

The core of the above implementation of the Deutsch-Jozsa 
algorithm is represented by a unitary, controlled black-box 
operator Uf applied between the Fourier transformed reg- 
ister and target CVs. Here, the Fourier transformed target 
CV together with a black-box operator induces a phase shift, 
which depends on the global property of the function f(x): 
Uf(\x)F t \7i/2)) = e- 2i f^\x)F t \7i/2) = e - in fW\x)F t \n/2). 
Notice that the state of the target CV is not changed follow- 
ing the action of Uf. In fact F t \n/2) is an eigenstate of Uf 
with an eigenvalue e~ lK f^ "kicked back" in front of the regis- 
ter C V [9] . Conventionally, the Deutsch-Jozsa algorithm em- 
ploys multiple quantum systems, however, as the above sim- 
ple analysis of the action of Uf indicates the target CV can 
be omitted. It is easy to show that a single register CV to- 
gether with a redefined black-box operator Uf = e~ 2lK l 2 f^ is 
enough to implement the Deutsch-Jozsa algorithm over con- 
tinuous variables. In Ref. [7], the authors arrived at the same 
conclusion, however, they used a slightly different approach. 
We emphasize that a direct consequence of employing a single 
system is that this protocol does not use any entanglement to 
determine the global property of the function in a single run. 
Moreover, the above implementation of the Deutsch-Jozsa al- 
gorithm is expressed in terms of the idealized position eigen- 
states. However, a more realistic and physically meaningful 
representation of a continuous variable is given by, for exam- 
ple, Gaussian states. 

Similar to the setting of discrete quantum systems (e.g. 
qubits), some features of the Deutsch-Jozsa algorithm can 
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FIG. 2: The general parameter estimation procedure involving state 
preparation P, evolution U((p) and generalized measurement M with 
outcomes x, which produces a probability distribution p(x\(p). 



FIG. 3: A quantum circuit representing the general protocol over 
continuous variables. The quantum network consists of the Fourier 
transform F and black-box gate Uf((p) applied to a single register 
CV prepared in the Gaussian state \G(xq)). The last operation is an 
inverse Fourier transformation that enables the interference of differ- 
ent computational paths. 



serve as a starting point for developing other quantum al- 
gorithms. A slightly modified black-box operator Uf = 
e -2i7i/2f(x) f or a s i m piifl e( i Deutsch-Jozsa algorithm can be 
used as the core of a protocol capable of estimating an un- 
known parameter that under appropriate conditions still re- 
tains the capabilities of the Deutsch-Jozsa algorithm. Before 
introducing this protocol, let us recall some basic concepts in 
quantum parameter estimation theory. 



III. PARAMETER ESTIMATION 

The most general parameter estimation procedure is shown in 
Fig. 2, and consists of three elementary steps: (i) prepare a 
probe system in an initial quantum state p(0), (ii) evolve it 
to a state p((p) by a unitary evolution U((p) = exp(— i(pJ4f), 
(iii) subject the probe system to a generalized measurement 
M, described by a Positive Operator Valued Measure (POVM) 
that consists of elements E x , where x denotes the measure- 
ment outcome. Here, the Hermitian operator is the gener- 
ator of translations in <p, the parameter we wish to estimate. 
The amount of information about (p that can be extracted by a 
measurement of the probe system is given by the Fisher infor- 
mation 



i 



p(x\q>) V dip 



dp(x\q>) 



(9) 



where p(x\(p) = Tv[E x p((p)] is the probability distribution 
given by the Born rule that describes the measurement data, 
and x is a discrete measurement outcome. Based on the Fisher 
information one can bound a minimal value of the uncertainty 
in <p with the quantum Cramer-Rao bound [10-12] 



(<5<p) 2 > 



1 



TF(cp) ' 



(10) 



where ( 5 (p) 2 is the mean squared error in the parameter (p, and 
T is the number of times the procedure is repeated. The ulti- 
mate limit of the quantum Cramer-Rao bound depends on how 
the Fisher information is bounded from above. The Fisher in- 
formation can be bounded in two ways: by the variance of J4? 
[13] or by the expectation value of J(? [14, 15] 



F(<p) < 16(A^) 2 and F(<p) < 4(J^) 2 



(11) 



where we again used h = \ . Since both bounds are completely 
general and complement each other, any parameter estimation 
procedure must respect them. Typically, the Fisher informa- 
tion is related to an appropriate resource count such as the 



average photon number, the average energy of the probe sys- 
tem or the number of unitary evolution gates that are used in 
the estimation procedure. The expectation value of ffl plays 
the role of the resource count [14]. We usually consider two 
scaling regimes of the quantum Cramer-Rao bound. The first 
regime: the so-called standard quantum limit (SQL) [16] or 
shot-noise limit is obtained when the Fisher information is a 
constant with respect to T and the resource count. The SQL is 
typically given by 



S(p 



> 



1 



(12) 



The second regime: the so-called Heisenberg limit [17] is ob- 
tained in a single- shot experiment (7 = 1) when the Fisher 
information scales quadratically with the resource count. The 
Heisenberg limit is then given by 



S(p> 



1 



(13) 



Therefore, the uncertainty in the parameter (p scales linearly 
inversely with the resource count. Both scaling regimes of 
the quantum Cramer-Rao bound can be compared directly in 
terms of an appropriate resource count [14]. 



IV. GENERAL PROCEDURE WITH GAUSSIAN STATES 

In this section, we present a general procedure capable of de- 
termining the value of a single parameter (p G [0, lit) or im- 
plementing the Deutsch-Jozsa algorithm (see Fig. 3). Here, 
the black-box operator is defined in the following way 



tf/(p) = exp(-2ip/(f)), 



(14) 



where f(x)\x) = f(x)\x). The function f(x) again takes only 
two values and 1. Without loss of generality, ideal, nonnor- 
malizable continuous variable states are regularized to Gaus- 
sian input states. Similar to the case of the Deutsch-Jozsa 
algorithm, any physical continuous-variables parameter esti- 
mation protocol can be implemented only on a finite domain. 
Therefore, we introduce the semi-Gaussian input state defined 
on a finite domain given by 



|G(*o)> 



f T dx 



(x-xq) 2 
2A 2 



(15) 



where A is the variance of the state and is the normalization 



constant given by N% 



erf(^a)+erf(^a 
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We note that for A <C T w e rec over the normalization con- 
stant in the form of N 2 = V 71 A 2 which is characteristic for a 
Gaussian state defined on an infinite domain, i.e. from — oo to 
+00. The Fourier transformed semi-Gaussian state defined on 
a finite domain can be written as 



\G(po)) = T XT ex P [-^(P-Po) 2 } \p), (16) 
J-pN p 



where 1 / (2A) is the variance of the Fourier transformed semi- 
Gaussian state and N p is given by 



Ni 



^n/AA 2 



[erf (2(P + po)A) + erf (2(P — po)A)] . 



For P ^> 1 / (2A) the normalization constant takes the form 
of Np = ^n/4A 2 , characteristic for a Fourier transformed 
Gaussian state define on an infinite domain. The relation- 
ship between domains of the semi-Gaussian input state and 
its Fourier transformed counterpart is given by P = 1/ (2T). 

The general procedure consists of the following instruc- 
tions: (i) prepare the register CV in the normalized semi- 
Gaussian state |r) = \G(xo)), and apply the Fourier transform 
F defined by 

F\x) = \x) p = -±=J T _ T dye 2i xy\y), (17) 

where \x) p is the Fourier transformed position eigenstate, i.e. 
the momentum eigenstate; (ii) subsequently, a black-box op- 
erator Uf((p) is applied. Then the state of the system is 



f 1 dx 

U f (q>)F\r) = -exp 
J-T N x 

= —f 



(x-xp) 
2A 2 



-2i<pf(x) I 



dxdy 



exp 



(X-Xp) 

2A 2 



xe' 



2iyx-2i(pf(y) 



\yh 



(iii) finally, an inverse Fourier transform F~ l is applied fol- 
lowed by a measurement. The state of the register CV is mea- 
sured by projecting onto the original semi-Gaussian state cen- 
tered around jcq. Measurement is described by a POVM set 
{Px ,Px }, where 



Px = J dxdyg xy \x) (y\ , and P XQ = I - P XQ 



with 



8 ^ = Nl 6XP 



(x-Xq) 
2£ 2 



exp 



(y-xof 

2e 2 



(18) 



(19) 



and £ is the intrinsic precision of the measurement apparatus, 
i.e. any continuous-variable measurement must have finite 
precision if it is to be physical, and N e is the normalization 

constant given by N 2 = \fn£ 2 j2 [erf ( ^ ) + erf ( ^ ) . 

The optimal measurement which corresponds to the initial 
semi-Gaussian register state has £ = A, thus N £ = N x . 



Now let us calculate the measurement statistics. Analyti- 
cal expressions for the measurement statistics are hard to find 
due to the presence of error functions erf(x). However, for 
the semi-Gaussian states with A <C T the calculations simplify 
considerably. Under this regime, the limits of integration for 
the integrals containing terms that depend on A range from 
— oo to +°°. Necessarily, the normalization constants have to 



be changed and are expressed as V2TN X = \f% v nA 2 . In other 
words, a semi-Gaussian input state defined on a finite domain 
is approximated with a Gaussian state defined on an infinite 
domain. Therefore, the measurement statistics based on the 
above POVM is given by the following expression 



p(xo\(p) 
p(x \(p) = 1 



4A- 



A z ^ 

71 j-p 

-p(x \(p). 



dzdye- 4A2 ^^e 2i ^-^ 



(20) 



Here, the interval (—P,P) is a finite domain of the Fourier 
transformed semi-Gaussian state \G(xo)) 9 and denotes the in- 
terval, where for this particular procedure function f(x) is de- 
fined. 

At this point, we have to give an explicit definition of the 
function. Functions f(x) defined on a finite domain return- 
ing only two values ({0, 1}) fall into three distinct categories: 
constant, balanced and neither constant nor balanced. We re- 
call that the objective of the Deutsch-Jozsa algorithm is to 
probe whether an unknown function f(x) is constant or bal- 
anced. We parameterize the three possibilities for defining 
f(x) by introducing a parameter r. The above integrals can 
then be evaluated for any function f(x) behaving as a step 
function, with the parameter r marking the point where f(x) 
changes its value. Hence, for r — and r = ±P the function 
f(x) is balanced and constant, respectively. For < r < P 
(or — P < r < 0), the function f(x) is neither constant nor 
balanced. We consider only positive values of r due to the 
symmetry of the setup. This leads to 

p(xo\q>) = ^[erf 2 (2PA)+erf 2 (2rA)] + 

^ [erf 2 (2PA) -erf 2 (2rA)] cos(2p), 
p(x \(p) = l-p(x \<p), 

where p(xo\(p) is the probability of measurement outcome to 
be in the interval xo ± £ and p(xo \ <p) is the probability of mea- 
surement outcome not to be in the interval xq ± £. 



A. Representations of f(x) 

Our choice to represent f(x) as a step function simplified our 
calculations. However, we can imagine more elaborate be- 
havior patterns for f(x). In principle, since in the case of the 
Fourier transformed idealized continuous variables all terms 
have amplitudes of equal magnitude, all finite sub-intervals, 
where the function takes value can be added up to a single 
interval. The same applies to all sub-intervals, where function 
takes value 1 . Therefore, one ends up with two intervals and a 
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relationship between them given by the parameter r. However, 
in the setting of semi-Gaussian states defined on a finite do- 
main, the above reasoning is not quite as straightforward. The 
amplitudes of the Fourier transformed Gaussian states have 
a slightly different magnitude. One may notice this feature 
by inspecting Eq. (20). Since in our calculations we favor a 
step-function representation over any other, let us estimate the 
maximum error we make with this assumption. Due to a triv- 
ial nature of a constant function, in the following analysis we 
consider a balanced function. We consider the step-function 
representation of a balanced function with r = 0. The biggest 
deviation from this representation is offered by a balanced 
function that changes its value twice at points r\ = —P/2 and 
V2=P/2. Both representations produce two distinct probabil- 
ity distributions PstepfaW) an d Phat{ x oW)i respectively, that 
differ by the error £p& given by 



£ PA = |1 -cos(2<p)| x 



24 



71 



(/>A) 6 + — (PAf + 0((PA) 

71 



10\ 



The error tends to zero with PA — » 0. This is natural since 
when A —> all amplitudes of the Fourier transformed ide- 
alized position eigenstate have the same magnitude, i.e., the 
spectrum is flat. 



B. Analysis 

Our procedure can be analyzed in two ways. As expected, 
from one perspective it behaves as a parameter estimation pro- 
tocol. From the other, it behaves as the Deutsch-Jozsa algo- 
rithm. First, we analyze the behavior of the parameter estima- 
tion part of the procedure. Based on the above measurement 
statistics, we calculate the Fisher information F((p). The min- 
imal value of F{(p) = occurs when function f(x) is constant 
(r = P) with the corresponding measurement statistics 



p(*o\<P) 
p(xo\<P) 



erf 2 (2/>A), 
l-erf 2 (2/>A). 



Conversely, the maximal value of the Fisher information 
4erf 2 (2/>A)(cos(2<p)-l) 



F(q>) 



erf 2 (2PA)(cos(2(p) + l)-2 



(21) 



occurs when function f(x) is balanced (r = 0) with the corre- 
sponding measurement statistics 

p(xo\<p) = ^erf 2 (2/>A)[l+cos(2(p)], 



M*o|<P) 



-erf 2 (2/>A)[l+cos(2(p)]. 



Here, the optimal value of the Fisher information F(<p) =4 
is given for erf 2 (2PA) = 1 P > 3/(2A) which, in general, 
implies P > 1/ (2A) and is consistent with the approximation 
applied above. The general dependence of the Fisher informa- 
tion F((p) on parameter r with P = 3/ (2A) and A = 1/ y/2 (the 
variance of the coherent state) is shown in Fig. 4. The dips that 
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FIG. 4: General dependence of the Fisher information F(q>) for five 
values of the parameter r: r = corresponds to the uppermost solid 
line (green), r = P/S corresponds to the dashed line (blue), r = P/4 
corresponds to the dashed-dotted line (brown), r = P/2 corresponds 
to the long dashed line (gray), and r = P corresponds to the lower- 
most solid line (red). 



are especially visible for the balanced function appear because 
the Fisher information F(<p) retains some dependence on the 
parameter <p since for P = 3/(2A): erf 2 (2PA) « 1. Based on 
the general dependence of F(<p) on r, we conclude that the 
maximal value of the Fisher information is indeed obtained 
for a balanced function. 

To address the optimality of our parameter estimation pro- 
tocol, we analyze the behavior of the generator of translations 
in the parameter <p: = f(x). The expectation value of the 
generator Jt? in the state of the register CV preceding applica- 
tion of the black-box operator, i.e. | Yin) = F\r) with A«T, 
is given by 

<JT> = </(*)> = l(erf(2PA) - erf(2rA)). 

Since f 2 (x) = f(x) the variance of the generator Jtf in | !//;•„) 
can be written as 



(Ajtry = (Am) 



^(erf(2PA)-erf(2rA)) x 
l-^(erf(2PA)-erf(2rA)) 



The maximal expectation value of the generator Jt? occurs for 
a balanced function (r = 0) with P > 3/ (2A) and is given by 
(Jtf) = 1/2. On the other hand, the maximal variance of the 
generator Jj? is (AJ^) 2 = 1/4. Hence, the Fisher information 
is bounded by F((p) < 16(AJ4?) 2 = 4. Therefore, we note 
that according to Eqs. (11) and (13) our procedure attains the 
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scaling regime of the Heisenberg limit. However to establish 
its optimality we must calculate whether S(p = 1/ y/F(q>). We 
use the standard expression for the mean squared error given 
by 

5 *'wmv <22) 

where X is the measurement observable defined as X = 
P XQ [see Eq. (18)]. Hence, for the final state W<p) = 
F~ l U f ((p)F\r) with£ = Awe calculate (X) = (VfylftolVfy) = 
\ erf 2 (2PA) [1 + cos(2<p)] . Based on the property P 2 Q = P XQ we 
find that (X 2 ) = (X). For P > 3/(2A) the mean squared error 
is 5 (p = 1/2. Hence, we conclude that for a balanced function 
our parameter estimation procedure over continuous variables 
attains the ultimate limit of the quantum Cramer-Rao bound, 
and therefore is optimal. This result constitutes an analogy to 
the phase estimation with a qubit realized as a single photon 
placed in the arms of the Mach-Zender interferometer. Here, 
the balanced property of function f(x) plays a role of two dis- 
tinct paths in a balanced Mach-Zender interferometer. 

Next, let us analyze the Deutsch-Jozsa side of the proce- 
dure. Under appropriate conditions the developed procedure 
can determine the character of function f(x). If a value of the 
parameter <p is fixed: (p = 7l/2 then the measurement statistics 
is given by 

p(xo) = erf 2 (2rA), 
p(xo) = l-erf 2 (2rA), 

It is clear that for a constant and balanced function f(x) the 
corresponding measurement statistics of the Deutsch-Jozsa al- 
gorithm are recovered. Indeed, when function f(x) is constant 
(r P) then 

p(xo) = erf 2 (2/>A), 
p(xo) = l-erf 2 (2/>A), 

and when function f(x) is balanced (r = 0) then p(xq) = 
and p(xo) — 1- The Deutsch-Jozsa algorithm over the semi- 
Gaussian states defined on a finite domain becomes a prob- 
abilistic procedure. This is consistent with the conclusions 
found in Ref. [7]. However, when the size of the domain is 
sufficiently large with P > 3/ (2A) then a definite distinction 
between constant and balanced functions can be made. Never- 
theless, even for large enough domains this implementation of 
the Deutsch-Jozsa protocol does not offer an unphysical, infi- 
nite speed-up over the classical procedures. We note that for 
ideal, nonnormalizable position eigenstates (A — > 0), the con- 
stant function measurement statistics is retained for P — > oo 
rendering P and r unphysical, thus making a meaningful dis- 
tinction between the balanced and constant functions impos- 
sible. 

We also calculated the Fisher information F(r) and plot- 
ted it against r G (0,P) for five different values of the param- 
eter (p = {^/2,57r/12,^/3,^/4,7r/8} with P = 3/(2A) and 
A = 1 / y/2 (see Fig. 5). The maximal value of the Fisher infor- 
mation F{r) is obtained for (p = 71 / 2 corresponding to a sim- 
plified Deutsch-Jozsa algorithm. We note that the optimality 
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FIG. 5: General dependence of the Fisher information F(r) for four 
values of the parameter cp: (p = n/2 corresponds to the uppermost 
solid line (green), (p = 5n/\2 corresponds to the dashed line (blue), 
cp = 7r/3 corresponds to the dashed-dotted line (brown), cp = n/4 
corresponds to the long dashed line (gray), and <p = 7t/8 corresponds 
to the lowermost solid line (red). The optimal value of r shifts from 
balanced to constant. 

changes from balanced to more constant when (p ^ n/2. Any 
further analysis of this side of the procedure is problematic 
due to a lack of the generator of translations in r. 

One possible application of the Deutsch-Jozsa part of our 
procedure is to test the quality of the implementation of func- 
tion f(x) employed in the parameter estimation protocol. 
Whenever the function is balanced or constant the quality of 
its implementation can be established by probing the parame- 
ter r. We also would like to stress that since we are employ- 
ing a single continuous variable, no entanglement is present 
at the preparation stage nor is created during the computation. 
The quantum superposition principle itself is responsible for 
speed-up over any classical procedure. Even though, in princi- 
ple, a single continuous variable is quite sufficient, a practical 
implementation of the Deutsch-Jozsa algorithm may require 
more continuous variables. 



V. CONCLUSIONS 

In conclusion, we developed a general procedure capable of 
performing two distinct tasks. For one mode of operation 
the protocol estimates a value of an unknown parameter with 
Heisenberg-limited precision. On the other hand, for a fixed 
value of the parameter in question the procedure addresses the 
Deutsch-Jozsa problem in a single run. Our procedure em- 
ploys Fourier transforms and black-box unitary operator ap- 
plied to a single continuous variable represented as the semi- 
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Gaussian state defined on a finite domain. Consequently, for 
this setup, the parameter estimation side of the procedure 
is optimal and the Deutsch-Jozsa algorithm offers finite, i.e. 
physically feasible, speed-up over any classical procedure. 
Furthermore, no entanglement is present at any stage of the 
procedure. A similar conclusions concerning quantum metrol- 
ogy can be found in Refs. [18, 19]. We emphasize a special 
role played by the balanced kind of functions f(x). The pro- 
cedure equipped with the black-box operator that introduces 
the parameter (p via the balanced function attains the ultimate 



limit of the quantum Cramer-Rao bound. This behavior can be 
linked to the phase estimation with a qubit realized as a single 
photon placed in the arms of the Mach-Zender interferometer. 
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